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ABSTRACT 



We describe a high-precision Doppler search for giant planets orbiting a well- 
defined sample of metal-poor dwarfs in the field. This experiment constitutes a 
fundamental test of theoretical predictions which will help discriminate between 
proposed giant planet formation and migration models. We present here details 
on the survey as well as an overall assessment of the quality of our measurements, 
making use of the results for the stars that show no significant velocity variation. 
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Subject headings: planetary systems: formation — stars: statistics — stars: 
abundances — techniques: radial velocities 

1. Introduction 

With a present-day catalogue of ~ 180 extrasolar planets 1 , several important statistical 
properties of the sample are beginning to emerge (for a review, see for example Marcy et al. 
2005 and Udry et al. 2006). Some of the most intriguing features unveiled so far include gas 
giant planets in few-day orbits, eccentricities spanning the entire range of possible values, 
the existence of a 'Brown Dwarf Desert', correlations among planetary orbital and physical 
parameters, evidence for a wide variety of dynamical interactions in multiple-planet systems, 
and a dependence of the frequency and properties of planetary systems on some of the 
characteristics of the parent stars (mass, metallicity). 

The largely unexpected properties of extrasolar planets have boosted theoretical research 
in the field of planetary system formation and evolution. Recently, improved understanding 
of some of the complex aspects of the underlying physical processes governing planet forma- 
tion in gaseous protoplanetary disks has enabled us to move from attempts to explain already 
known features of the exoplanet sample to more refined models making testable predictions. 
In turn, surveys for planets with a variety of techniques have started targeting more carefully 
defined and selected stellar samples, with the aim of providing the observational evidence 
required to address a specific scientific question (rather than simply focusing on the discovery 
of extrasolar planets), and thus helping to discriminate between different theoretical planet 
formation models. 

This is the first of a series of papers in which we present results from a spectroscopic 
search for giant planets orbiting a sample of metal-poor stars in the field. The two competing 
mechanisms for gas giant planet formation, core accretion and disk instability, produce sig- 
nificantly different distributions of planet masses and orbital elements (e.g., Rice et al. 2003; 
Ida & Lin 2004a), and they predict a very different dependence of planetary frequency on 
stellar metallicity and mass (Boss 2002, 2006; Ida & Lin 2004b, 2005; Laughlin et al. 2004). 
In order to determine whether core accretion or disk instability is the dominant formation 
mode for giant planets, or to verify the existence of bimodal planet formation, one should 
then compare the frequency of gas giant planets and their properties between metal-rich 
and metal-poor stars. However, the low-metallicity stellar sample that has been searched 
for planets is at present too small to test but the most outstanding differences between such 
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hypothetical populations. It is now crucial to provide a statistically significant, unbiased 
sample of metal-deficient stars screened for giant planets. 

We have arranged this first paper as follows. In Section 2 we present the scientific 
case for a Doppler survey for giant planets orbiting field metal-poor dwarfs. We describe 
in Section 3 some of the technical aspects of our Keck/HIRES radial-velocity survey. We 
present in Section 4 preliminary results from the first 3 years of observations, focusing on the 
detailed assessment of the quality of our measurements, making use of the results for all the 
stars that show no significant velocity variation. We provide in Section 5 a brief summary 
and conclusions. 



2. Testing Giant Planet Formation and Migration Models 

2.1. Stellar Metallicity and Planets: Observations 

The connection between the presence of giant planets and the metal content of the 
parent stars has been the subject of a significant number of studies in the past (see for 
example Gonzalez 2003, and references therein). The average of the metallicity distribution 
of planet-hosting stars is [Fe/H] ~ 0.14 (e.g., Fischer & Valenti 2005), whereas the mean 
value for the solar neighborhood is [Fe/H] ~ —0.1 (Nordstrom et al. 2004). With improved 
statistics, in recent years the hypothesis that super-solar metallicity could correspond to a 
higher likelihood of a given star to harbor a planet has been conclusively proved (e.g., Santos 
et al. 2001, 2004; Fischer & Valenti 2005). As of today, the frequency of giant planets around 
metal-rich ([Fe/H] > 0.3) dwarfs (integrating over the F-G-K spectral types) is f p > 20%, 
while this fraction decreases to f p ~ 3% for metal-poor stars (—0.5 < [Fe/H] < 0.0). For 
metallicities below solar, f p appears to be roughly constant (Santos et al. 2004a; Fischer & 
Valenti 2005). 

Many authors have debated whether the observational evidence is an indicator of pri- 
mordial high metallicity in the planet host stellar sample (Santos et al. 2001; Reid 2002), or 
if the trend with [Fe/H] could be due to pollution by ingested planetary material (Laughlin 
2000; Gonzalez et al. 2001; Israelian et al. 2001; Murray & Chaboyer 2002). The idea of 
pollution is losing credit among the scientific community, primarily based on the lack of 
any correlation between [Fe/H] and effective temperature T c s, or convective envelope mass 
M conv , for the planet host sample (e.g. Pinsonneault et al. 2001; Fischer & Valenti 2005; 
Santos et al. 2003, 2004a, and references therein). Results on this specific issue are not yet 
conclusive, however. For example, Vauclair ( 2004) has recently pointed out how the absence 
of a [Fe/H]-M conv correlation does not automatically imply that stars with planets have not 
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been polluted. 

Furthermore, theoretical calculations (Montalban & Rebolo 2002; Boesgaard & King 2002) 
suggest that detection of anomalous abundances of rare elements such as lithium (Li) or 
beryllium (Be) could be interpreted as evidence for recent accretion of planets onto the at- 
mosphere of a star. The abundances of Li isotopes in the spectral region around the 6707A 
line in planet-host stars have been investigated in the recent past by several authors (Gonza- 
lez & Laws 2000; Ryan 2000; Israelian et al. 2001, 2003, 2004; Reddy et al. 2002; Mandell et 
al. 2004), and similar studies have been conducted for the Be II lines at 3130A and 313lA. 
(Garcia Lopez & Perez de Taoro 1998; Deliyannis et al. 2000; Santos et al. 2002, 2004b). 
While the presence of the 6 Li isotope has actually been detected in some planet-harboring 
stars (Israelian et al. 2001, 2003; Laws & Gonzalez 2001), suggesting that accretion of plane- 
tary material can indeed take place in some stars, in general stars with planets have normal 
light-element abundances, typical of field stars. It thus seems unlikely that pollution effects 
can be responsible for the overall metallicity enhancement of the planet host stellar sample. 

In addition, analyses of over a dozen other elements have been carried out in the recent 
past (Santos et al. 2000; Gonzalez et al. 2001; Smith et al. 2001; Takeda et al. 2001; Sadakane 
et al. 2002; Zhao et al. 2002; Bodaghee et al. 2003; Ecuvillon et al. 2004a, 2004b, 2006a; 
Beirao et al. 2005; Gilli et al. 2006; Sozzetti et al. 2006), and the general evidence is that the 
abundance distributions in stars with planets are the extension of the observed behavior for 
[Fe/H], a result quantified by trends of decreasing [X/Fe] with decreasing [Fe/H]. The absence 
of any statistically significant trend of metallicity [X/H] with condensation temperature T c 
(e.g., Sozzetti et al. 2006, Ecuvillon et al. 2006b; Gonzalez 2006) is one more piece of 
circumstantial evidence that the best explanation for the metallicity excess in stars with 
planets is that the enhanced [Fe/H] is primordial in nature. 

Several studies have also focused on possible correlations between stellar metallicity and 
planet properties. While no significant trend was found between [Fe/H] and planet mass or 
orbital eccentricity (e.g., Udry et al. 2002; Santos et al. 2001, 2003; Fischer et al. 2002), 
evidence (albeit weak) for a correlation between the metallicity of planet-harboring stars and 
the orbital periods P of the planets has been pointed out (e.g., Sozzetti 2004, and references 
therein; Santos et al. 2006). This correlation is highlighted by an excess of close-in planets, 
on few-days orbits, around the metal-rich ([Fe/H]> 0.0) sample of planet hosts. 
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2.2. Stellar Metallicity and Planets: Theory 

Within the framework of the scenario of gas giant planet formation by core accretion 
(e.g., Lissauer 1993; Pollack et al. 1996; Alibert et al. 2005), recent studies have successfully 
reproduced the strong dependence of planetary frequency on stellar metallicity, in qualita- 
tively good agreement with the observed trend (Kornet et al. 2005; Ida & Lin 2004b). The 
probability of forming giant planets according to the disk instability model (e.g., Boss 1997, 
2000; Mayer et al. 2004), however, is remarkably insensitive to the primordial surface density 
of solids of the protoplanetary disk (Boss 2002; Rice et al. 2003), thus planet occurrence 
should not be hampered around metal-poor stars. On the one hand, the observed trend sug- 
gests that giant planet formation by core accretion predominates in the metal-rich regime 
([Fe/H]> 0.0). On the other hand, f p appears to be rather flat in the metal-poor regime 
([Fe/H]< 0.0). The possible evidence for bi-modality of the planet frequency distribution as 
a function of metallicity (Santos et al. 2004; Fischer & Valenti 2005) suggests the existence 
of two different mechanisms for forming gas giant planets. However, due to the low numbers 
of metal-poor stars ([Fe/H]< —0.5) surveyed to-date, no definitive conclusion can be drawn, 
except that maybe both mechanisms operate (Beer et al. 2004). 

If real, several possible explanations can be put forth for the existence of a P— [Fe/H] 
correlation. For example, migration rates might be slowed down in metal-poor protoplane- 
tary disks (Livio & Pringle 2003; Boss 2005), although the predicted dependence of migration 
time-scales on [Fe/H] appears to be somewhat weak. The correlation may also arise as a 
consequence of longer time-scales for giant planet formation around metal-poor stars, and 
thus reduced migration efficiency before the disk dissipates (Ida & Lin 2004a; Alibert et al. 
2005). Another possibility (Santos et al. 2006) is related to planetary internal structure 
arguments: if planets formed in low-metallicity disks have small rocky cores (e.g., Pollack 
et al. 1996; Ida & Lin 2004), their low density might hamper survival against evaporation 
(Baraffe et al. 2004; Lecavelier des Etangs et al. 2004), once they have migrated to very 
close-in orbits. However, even in this case small-number statistics in the low-metallicity 
regime prevents one from reaching a clear conclusion. 

One way or another, in order to unambiguously determine the role of metallicity in 
gas giant planet formation (for example, is f p a truly monotonic function of [Fe/H], or is 
planetary frequency constant in the metal-poor regime?), and consequently discriminate 
between proposed explanations for the observed trends in the data, it is crucial to provide a 
statistically significant, unbiased sample of metal-poor stars screened for giant planets. 
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2.3. Searching for Planets Around Field Metal-Poor Stars 

The absence of short-period (P < 8.3 days) transiting planets in the moderately metal- 
deficient ([Fe/H] ~ —0.7) globular cluster 47 Tucana? has been used by Gilliland et al. (2000) 
and by Weldrake et al. (2005) to argue that low-metallicity stars are less likely to harbor 
giant planets. The claims by these authors suffer however from some ambiguity, because 
in the cluster core investigated by Gilliland et al. (2000) with HST transit photometry, 
crowding could play a significant role in giant planet formation, migration, and survival 
(e.g., Davies & Sigurdsson 2001; Bonnell et al. 2001; Hurley & Shara 2002; Fregeau et al. 
2006). The outer regions of the cluster monitored by Weldrake et al. (2005) are less affected 
by crowding. However, the lower occurrence rate of hot Jupiters in a metal-poor environment 
does not rule out the existence of a population of giant planets at wider radii. Indeed, this 
possibility has recently been supported by the findings of Sigurdsson et al. (2003), who, 
using HST data, were able to infer a mass of a few Jupiter masses for the third, long-period 
component orbiting the white dwarf - pulsar system B1620-26 in the globular cluster M4, 
five times more metal-poor than 47 Tuc. Their results provide the first evidence for planet 
formation in very metal-poor environments. In light of the M4 announcement, the 47 Tuc 
results could be re-interpreted as follows: Dynamical disruption in dense clusters is not 
sufficient to completely destroy any planetary population, and the lack of transiting planets 
on short-period orbits might be due to other processes, such as a metallicity dependence 
in the migration mechanism, or suppression of migration (but not formation) in globular 
clusters 2 . 

By addressing the field population of metal-poor stars directly, it is then possible to 
eliminate dynamical interactions in dense stellar environments as a possible source of inter- 
ference with planet formation, or with migration to close-in orbits, or with planet survival. 



2 Notably, the lack of any apparent correlation in the plane defined by minimum mass M p sini and orbital 
period P for planets in binary and triple stellar systems (as opposed to the observed M p sin i — P correlation 
in the case of planets orbiting single stars) is interpreted as evidence for enhanced migration efficiency for 
planets formed in stellar systems (Zucker & Mazeh 2002; Udry et al. 2003; Eggenberger et al. 2004). 
However, the birthplaces for such systems are stellar groups and clusters with much lower stellar densities 
than globular clusters. Thus, both the impact on planet formation and migration efficiency as well as the 
evolutionary history of dynamical interactions are likely not the same in these two different environments 
(e.g., Davies & Sigurdsson 2001; Adams et al. 2006) 
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3. The Keck/HIRES Doppler Survey of Metal-Poor Dwarfs 

In 2003, we began a high-precision radial-velocity survey of ~ 200 metal-poor stars, 
using HIRES on the Keck 1 telescope (Vogt et al. 1994) and its I2 gas absorption cell as the 
reference velocity metric (Butler et al. 1996). The goal of this project is to carry out an 
initial reconnaissance for gas giant planets orbiting within 1 AU of a statistically significant 
sample of low-metallicity dwarfs. 

The sample has been drawn from the Carney-Latham and Ryan samples of metal- 
poor, high-velocity field stars (e.g., Carney et al. 1994; Ryan 1989; Ryan & Norris 1991). 
The initial combined database totaled 1558 objects. A number of selection criteria have 
been adopted in order to finalize our list of targets. First of all, a key advantage of the 
Carney-Latham and Ryan samples is that we have monitored the radial velocities of the 
stars in these samples, mostly for more than 3,000 days, using the CfA Digital Speedometers 
(Latham 1992). This has allowed us to identify most of the stars with stellar companions 
that would interfere with the formation or survival of planets in the habitable zones (Carney 
et al. 2001; Latham et al. 2002). All stars included in the final list of targets showed no 
sign of velocity variation at the 0.5 to 1.0 km s~ x level. Second, significant chromospheric 
activity (quantified for example through the chromospheric emission ratio -Rhk) anc ^ l &r g e 
values of stellar rotational velocity v sini should be avoided, as they constitute sources of 
intrinsic radial velocity "jitter" that can mask, and sometimes even mimic, the presence of 
orbital reflex motion due to planetary mass companions (Saar et al. 1998; Santos et al. 2000; 
Queloz et al. 2001; Paulson et al. 2004). Fortunately, old stars have the advantage that they 
rotate slowly and have low levels of chromospheric activity. All of the stars in our sample 
exhibit rotational velocities v shu < 10 km s _1 , and most of them have rotations below the 
resolution limit of the CfA Digital Speedometers (8.5 km s _1 ), so that a value of v sini could 
not be determined. Thus, we do not expect velocity jitter due to astrophysical phenomena 
to be a problem for this sample. 

However, metal-poor stars have weak absorption lines in comparison to their solar- 
metallicity counterparts. The lines also grow weaker as the effective temperature rises. 
Furthermore, very metal-poor stars are rare, and therefore they tend to be distant and 
faint. In order to characterize the behavior of the radial velocity precision orv as a function 
of stellar metallicity, effective temperature T eff , and visual magnitude V (assuming non- 
rotating, inactive stars), we have run simulations utilizing the CfA library of synthetic stellar 
spectra (e.g., Nordstrom et al. 1994; Latham et al. 2002). We show in Figure 1 how <trv 
degrades as a function of [Fe/H] and T e ff, assuming a fixed exposure time and a typical 
measurement precision of 10 m s _1 for [Fe/H] = 0.0. For a solar-type star, [Fe/H] = —1.0 
corresponds to a degradation in ory of a factor ~ 2, while for a significantly cooler star, with 
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a more complex spectrum 3 , the effect is less severe (in addition, surface gravity constitutes 
only a higher-order effect). This result confirms the empirical findings of Santos et al. (2003) 
and Fischer & Valenti (2005), based on the analyses of the achieved velocity precision as 
a function of metallicity carried out with the stellar databases of their respective Doppler 
planet surveys. Those studies were aimed at ruling out possible observational biases that 
might contribute to the observed correlation between f p and [Fe/H]. In both cases, Santos et 
al. (2003) and Fischer & Valenti (2005) conclude that, given the typical single-measurement 
precision a RV ~ 3 — 5 m s -1 achieved on bright objects (typically V < 8.0), a velocity 
degradation of a factor 1.5 — 2 does not imply that a fraction of the giant planets orbiting low- 
metallicity stars might have gone undetected (at least in the orbital period range presently 
covered by Doppler surveys). However, not many bright metal-poor stars (—0.5 < [Fe/H] 
< 0.0) are found in the solar neighborhood (d < 50 pc). In order to create a statistically 
significant database (hundreds of stars) of metal-poor ([Fe/H] < —0.5) dwarfs, one must 
then include more distant, and fainter, objects. For the purpose of optimizing the trade- 
off between number of objects surveyed and total observing time required, we have refined 
our sample of metal-poor dwarfs from the Carney-Latham and Ryan surveys by adopting 
magnitude and temperature cut-offs (V < 12.0 and T cS < 6000 K, respectively), and by 
selecting objects in the metallicity range —2.0 < [Fe/H] < —0.6. 

In the empirical error model we have obtained using the simulations with the CfA 
database of synthetic spectra, the radial-velocity uncertainty is 



where F([Fe/H], T eS ) is an empirical function of temperature and metallicity based on the 
simulation results shown in Figure 1, and to and Vq are reference integration time and mag- 
nitude for a star with the temperature and metallicity of the Sun. Based on our experience 
with solar neighborhood G dwarfs observed with HIRES for the G Dwarf Planet Search 
Program (Latham 2000), we have set a threshold of 20 m s _1 precision for planet detection, 
and computed the relative exposure times needed to achieve such precision, for each star in 
our sample. Furthermore, we decided to limit the maximum exposure times to 15 minutes, 
to minimize uncertainties in the barycentric correction. We show in Figure 2 the distri- 
bution of [Fe/H], V, and distance estimates (for those objects with Hipparcos parallaxes, 
photometric otherwise) for the sample of 278 metal-poor stars derived after adopting all the 
selection criteria described above. In order to highlight the different ranges of metallicity, 



3 The velocity information content depends on the mean absolute value of the slope of the spectrum, which 
increases for later spectral types (e.g., Butler et al. 1996; Bouchy et al. 2001) 
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(1) 



- 9- 



magnitude, and distance spanned by our survey, the same distributions for a large fraction of 
the present-day sample of planet hosts are also shown. Given that the average stellar mass 
in our sample is ~ 0.69 M , setting a velocity precision threshold at 20 m s -1 is sufficient to 
detect (at the 5-cr level) velocity variations of planetary companions with minimum mass in 
the average range 0.59 Mj < M p sini < 2.75 Mj, or higher, for orbital periods in the range 
0.01 < P < 1 yr. 

Finally, if it turns out to be true that planets did not form around metal-poor dwarfs, 
then we need to observe a large-enough sample so that a null result, i.e. no detections, is 
significant. The rate at which giant planets with orbits inside 1 AU are being discovered by 
radial velocities appears to be about P ~ 3% (e.g., Marcy et al. 2005). In order for the 
failure to detect any planetary companions to be significant at the 3-<r level (corresponding 
to a probability of 0.0027), we need to survey a sample of N stars, where (1 — P) N = 0.0027, 
which is satisfied for N = 194. We therefore defined our final target list by selecting a subset 
of 200 metal-poor stars out of the abovementioned larger sample of 278 objects, which will 
eventually provide a robust 3-cr null result in the case of no detections. 



4. Preliminary Results 

Our observing program has been awarded an average of 2 Keck/HIRES nights per 
semester, starting in early 2003. Good temporal coverage is thus a serious issue, as such 
a scheduling of the observations is good for sampling long-period radial- velocity variations, 
but rather poor for identifying possible short-period variables. We have tried to obtain at 
least three velocity measurements per star per year, with at least one set of back-to-back 
observations (taken in two consecutive nights), in order to mitigate our bias toward poor 
sampling (and thus significant aliasing and ambiguities) at short orbital periods. 

The first important step is to provide an assessment of the long-term stability of the 
velocity zero-point and single-measurement precision achieved for planet detection, in light 
of the predicted exposure times needed to reach o^y = 20 m s" 1 . Our analysis pipeline 
incorporates the full modeling of temporal and spatial variations of the HIRES instrumental 
profile (Valenti et al. 1995), similarly to the method adopted by other groups (e.g., Butler 
et al. 1996; Korzennik et al. 2000; Cochran et al. 2002). Spectral modeling of each echelle 
order containing I 2 lines is carried out independently, and internal uncertainties for each 
observation are computed from the scatter of the velocities around the mean. This analysis 
technique has allowed us to significantly improve upon our initial estimates of achievable 
radial velocity precision. In Figure 3 we show the histogram of the rms velocity residuals 
for about 80% of our sample, for which we have achieved substantially uniform temporal 
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coverage (typically 5-6 observations per star, with at least two back-to-back exposures, over 
a time-span of at least a year). The rms velocity residuals distribution of the full sample 
(excluding variables with rms > 30 m s _1 ) averages ~ 9 m s" 1 . In Figure 4 we show the rms 
velocity residuals as a function of the time-span of the observations. Overplotted are the 
median and standard deviation in 500-days bins (again, variable stars are not included). As 
about two dozen of the stars in our sample are in common with the G dwarf planet survey 
of Latham (2000), we could establish the long-term stability of the velocity zero-point over 
time-scales of up to eight years. This demonstrates the true radial-velocity precision we are 
obtaining on the sample of metal-poor stars, with a significant improvement of over a factor 
of 2 with respect to the targeted 20 m s _1 single-measurement precision. 

The exposure times predicted by the model derived from the simulations with the CfA 
library of stellar spectra are determined as a function of [Fe/H], T cff , and V. However, our 
program stars are up to 3-4 magnitudes fainter and up to over 100 times more metal-poor 
than typical targets of Doppler planet searches. One possible matter of concern would then 
be the evidence of systematic trends in the velocity scatter as a function of these three 
parameters. However, as shown in Figure 5, no clear rms velocity trends as a function of 
[Fe/H], T e ff, and V are present. This gives us confidence that the model we developed for 
the dependence of the radial velocity precision on the above parameters is robust. 

A more quantitative, and challenging, test to demonstrate that the character of our 
errors is well understood can be carried out by studying the distribution of the velocity 
residuals compared with their formal uncertainties. We define the ratio ARV/ctrv as the 
difference between the velocity values and their mean value for each star, divided by their 
estimated uncertainties. In the ideal case, if internal errors are an accurate tracer of the 
actual uncertainties in the measurements, this ratio should have a Gaussian distribution 
with zero mean and unit dispersion. Discrepancies between predicted and actual errors 
should be reflected in measurable departures from gaussianity. We show in Figure 6 the 
histogram of the ratio ARV/ctrv for all our program stars (including variables). Overplotted 
is a reference Gaussian with zero mean and unit dispersion. The ARV/ctrv distribution is 
very close to Gaussian, with no apparent positive or negative biases. Slightly elevated wings 
(~ 6% of the velocity differences are larger than 5 sigma) indicate the presence of either 
non-Gaussian outliers or true variables. The result shown in Figure 6 demonstrates that 
our internal errors are realistic. The fundamental conclusion is that we achieve sufficient 
radial-velocity precision in our sample to clearly detect Jupiter-mass objects within 1 AU of 
metal-deficient dwarfs. 

We show in Figure 7 some examples of our velocity measurements of metal-deficient 
dwarfs. The top two panels show results for two of the stars (HD 157089 and G68-30) 
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in common between our survey and the G Dwarf Planet Search program (Latham 2000). 
These objects are constant to ~ 10 m s _1 over a time-scale of ~ 8 years. The examples 
of velocity time-series shown in the other panels of Figure 7 cover a range of ~ 1000 K in 
T e ff, ~ 4 magnitudes, and ~ 1 dex in [Fe/H]. The rms of the observations range between 6 
m s" 1 and 11 m s _1 , with an average of ~ 9 m s _1 . While a large fraction of our sample 
shows no significant velocity variations over the time-span of the observations, a number of 
objects do exhibit velocity variability indicative of the presence of companions. A thorough 
analysis of our planet detectability thresholds and a detailed presentation of all our velocity 
measurements will be the subject of our second paper. 

5. Conclusions 

One way to distinguish between proposed models of gas giant planet formation is to 
confirm or rule out, on an observational basis, their different predictions on planet frequency 
f p as a function of the metallicity of the stellar hosts. To address this issue, we have 
undertaken a Doppler search for giant planets within 1 AU of a sample of 200 metal-deficient 
(—2.0 < [Fe/H] < —0.6) dwarfs in the field. This is the sample size needed to provide a 
statistically significant result (at the 3 — o level) in the case of no detections. Using the 
Keck 1 telescope and its HIRES spectrograph, we have achieved a long-term radial-velocity 
precision of ~ 9 m s _1 , independently of [Fe/H] and T e g, and for stars 2 to 4 magnitudes 
fainter than the targets for most other radial-velocity planet surveys. We have provided 
convincing evidence that our internal errors estimates are reliable, and thus demonstrated 
that we achieve sufficient radial- velocity precision in our sample to clearly detect Jupiter-mass 
objects within 1 AU of metal-deficient dwarfs. A number of objects with significant radial- 
velocity trends have been identified, for which we plan to perform follow-up observations. 
At the conclusion of our Doppler survey, we will be able to place useful upper limits on the 
existence of planetary companions of given mass and period around metal-poor stars, and we 
will then compare the frequency of gas giant planets and their properties between metal-rich 
and metal-poor stars. These issues will be addressed in future papers. 
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Fig. 1. — Degradation in the radial velocity precision ctrv as a function of stellar metallicity, 
effective temperature, and gravity, for fixed exposure time. The zero point of ctrv for solar 
values of [Fe/H], T cff , and log 5 is arbitrarily scaled to 10 m s _1 . 

Fig. 2. — Distributions of visual magnitudes (top), distances from the Sun (center), and 
metallicities (bottom) for a sample of 278 metal-poor stars selected with the criteria detailed 
in the text. For comparison, the same distributions for a large sample of 119 planet hosts 
are also shown (data from Santos et al. 2004a, 2005, and Sozzetti et al. 2004). 

Fig. 3. — Rms velocity distribution for ~ 80% of the stars in our sample, for which uniform 
temporal coverage has been obtained (see text for details). Objects exhibiting significant 
radial velocity variations (> 30 m s^ 1 ) are not shown. 

Fig. 4. — Rms velocity residuals as a function of the time-span of the observations. Overplot- 
ted are the median (large filled circles) and standard deviation in 500-days bins. Variables 
with rms > 30 m s _1 are not taken into account. 

Fig. 5. — Radial velocity scatter (excluding variables with rms > 30 m s" 1 ) as a function of 
[Fe/H] (left), T c s (center), and V (right). Based on a rank-correlation test, the probability 
of no correlation in the three cases is 0.61, 0.14, and 0.09, respectively (i.e., no significant 
correlation is present). 

Fig. 6. — Histogram of all the velocity residuals (including variables) normalized by their 
formal uncertainties. The dotted line represents a reference Gaussian distribution with zero 
mean and unit dispersion. If, for example, internal errors are over- or under-estimated, this 
effect should show up as distortions in the distribution. The width of the distribution is very 
close to unity, indicating the absence of significant biases. Formal errors are a good estimate 
of the true underlying uncertainties. 

Fig. 7. — Observed relative velocities <5RV for a sample of stars in our program. The top 
two panels show objects with constant velocity to ~ 10 m s _1 over a time-span of 8 years, 
after combining observations taken in the context of the G Dwarf Planet Search program 
(Latham 2000). The mean rms of the stars shown in the other panels is 9 m s" 1 , the same as 
the average of the full sample (excluding variable stars). Comparable precision is achieved 
over ranges of 1000 K in T eff , 4 magnitudes, and 1 dex in metallicity. 
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